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At  the  time  of  formulating  the  statement  of  the  problem  to 


of  NAHBE  combustion  was 


as  defined  in  Ref.  1.  That  is.  a 


understanding  of  the  hypothesis  of  combustion  with  pressure 


exchange.  This  work  was  undertaken  with  the  hope  of  gaining 


some  quantitative  understanding  of  the  upper  and  lower  bounds 
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ABSTRACT 


This  report  describes  the  analytical  work  done  on  the 
Naval  Academy  Heat  Balanced  Engine  (NAHBE)  Combustion  Process. 
The  purpose  of  this  one  year  study  is  to  investigate  the  pro- 
posed utilization  of  pressure  exchange  mechanism  and  the  basic 
nonsteady  flow  processes  and  combustion  phenomena  associated 
with  the  NAHBE  combustion  system,  utilizing  the  experimental 
test  data  obtained  over  a period  of  years  at  the  U.S.  Naval 
Academy,  Annapolis.  A preliminary  analysis  of  the  effect  of 
pressure  exchange  on  the  performance  of  the  NAHBE  engine  is 
given  in  reference  (5)*  based  on  the  "volume  mode"  approach. 

A *front  mode''  approach  is  assumed  to  construct  a wave 
diagram  as  a first  step  towards  an  understanding  of  the  non- 
steady flow  dynamics  of  NAHBE  combustion  system.  Further,  two 
simple  heat  release  models  are  used  to  obtain  expressions  for 
pressure  and  the  results  compared  with  experimental  NAHBE 
test  data. 
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Specific  heat  ratio 


Constant 

Crank  angle  measured  from  top  dead 
center 

Equivalence  ratio 


Combustion 

Initial 

Shock 

Total 


I 

I 


INTRODUCTION 


1 


The  nature  of  Naval  Academy  Heat  Balanced  Engine  (NAHBE) 
Combustion  System  is  described  in  reference  1.  Because  of  the 
lack  of  information  about  the  details  of  the  actual  combustion 
process,  a preliminary  analysis  based  on  the  "volume  mode" 
approach  was  initiated  to  determine  the  effect  of  pressure 
exchange  on  the  performance  of  the  NAHBE  engine.  Fig.  1 shows 
the  basic  analytical  model  used  in  this  analysis.  Region  1 is 
occupied  by  the  combustion  gas,  region  2 by  the  balancing  gas. 
The  interface  between  the  two  is  simulated  by  a massless, 
frictionless,  impermeable,  and  adiabatic  piston.  Based  on  the 
volume  mode  analysis,  it  is  shown  that  the  p(t)  history  during 
combustion  is  independent  of  the  presence  or  absence  of  a 
balancing  chamber  both  for  the  case  of  stationary  power  piston 
and  the  more  general  case  of  a moving  piston.  Further,  the 
analysis  shows  that  for  both  cases,  the  proposed  utilization 
of  pressure  exchange  has  the  effect  of  increasing  the  available 
work  for  any  given  energy  input. 
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PRELIMINARY  NON-STEADY  FLOW  ANALYSIS  OF  NAHBE  COMBUSTION 
PROCESS  WITH  ADIABATIC  BALANCING  CHAMBER  AND 
STATIONARY  POWER  PISTON 

A "front"  mode  of  combustion  is  assumed  for  the  preliminary 
non- steady  flow  analysis,  wherein  it  is  assumed  that  the  chemi- 
cal reactions  take  place  instantaneously  as  the  unburned  gas 
passes  through  an  advancing  flame  front.  This  approach  con- 
siderably reduces  the  procedures  involved  in  the  construction 
of  the  wave  diagrams  on  either  side  of  the  flame  front.  The 
main  objective  of  the  construction  of  the  wave  diagram  is  to 
get  an  immediate  "feel"  for  the  nature  and  type  of  interactions 
and  transients  involved  and  their  relative  strengths.  Also 
one  can  approximately  compute  the  combustion  time  which  plays 
an  important  role  in  the  nature  and  formation  of  the  combustion 
products . 

Fig.  2 shows  the  analytical  model  used  for  the  non-steady 
flow  analysis.  As  shown  in  Figure  1,  region  1 is  occupied  by 
the  combustion  gas,  region  2 by  the  balancing  gas  (Air). 
Initially  the  interface  which  merely  separates  the  combustible 
mixture  of  gases  from  air  is  assumed  to  be  located  very  close 
to  the  balancing  chamber  at  the  beginning  of  ignition.  But 
as  long  as  the  interface  is  not  reached  by  the  flame  front, 
the  state  of  the  gas  is  the  same  on  both  sides. 

The  NAHBE  combustion  differs  from  that  of  a conventional 
spark  ignition  combustion  mainly  because  of  the  presence  of  a 
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secondary  chamber  called  the  balancing  chamber.  Reference  1 
points  out  the  fact  that  no  combustion  has  been  observed  in 
the  balancing  chamber  which  merely  acts  as  a reservoir  of 
hot  compressed  air.  According  to  the  pressure  exchange 
mechanism  proposed  therein, the  balancing  air  is  fed  into  the 
main  chamber  later  in  the  combustion  cycle  as  a result  of 
interactions  and  thereby  the  combustion  proceeds  to  completion 
over  a longer  period  of  time.  Therefore  the  movement  of  the 
interface  into  the  balancing  chamber  and  the  resulting  inter- 
actions are  of  great  importance  to  an  understanding  of  NAHBE 
combustion  process. 

The  plane  circular  interface  once  after  reaching  the 
balancing  chamber  moves  as  a cylindrical  interface  having  an 
area  given  by 


area  = 2ir  rt  (1) 


where 


r = instantaneous  radius 
t = balancing  chamber  length 

Note  that  the  area  of  this  cylindrical  interface  decreases 
with  time  as  the  interface  approaches  the  stem  of  the  balancing 
chamber  prior  to  being  reflected  back.  Hence,  the  movement  of 
the  interface  in  the  combustion  chamber  and  the  balancing 
chamber  can  be  visualized  as  taking  place  in  a conical  tube 
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with  a gradual  reduction  in  the  surface  area  of  the  interface. 

This  gradual  change  in  area  can  be  idealized  as  a finite 
number  of  discrete  change  in  area  for  the  purpose  of  construct- 
ing the  details  of  the  interactions  on  a wave  diagram.  This  is 
shown  in  Figure  3. 

Once  ignition  of  the  mixture  occurs  at  the  closed  end,  a 
flame  front  is  produced  and  advances  into  the  unburned  mixture. 
Since  the  flame  is  assumed  to  reach  its  burning  velocity  in- 
stantly upon  ignition,  the  acceleration  of  the  gas  takes  place 
through  a shock  wave  traveling  ahead  of  the  flame  front.  This 
results  in  a series  of  interactions  and  transients  which 
include : 

(1)  The  shock  wave  generated  at  the  closed  end  interacts  with 
the  first  discontinuous  change  of  cross  section. 

(2)  Deflection  of  the  interface  by  the  reflected  shock  wave 
at  the  first  change  of  cross  section. 

(3)  Interaction  of  the  transmitted  shock  wave  at  the  second 
change  of  cross  section,  etc. 

The  interactions  and  transients  that  subsequently  follow 
are  shown  in  the  wave  diagram  of  Figure  4. 

The  following  assumptions  are  made  in  the  non-steady  flow 
analysis : 

(a)  The  flame  front  advances  into  the  unburned  gas  with 
a burning  velocity  that  is  taken  proportional  to  the 
absolute  temperature  of  the  unburned  gas. 


(b)  Initial  air/fuel  ratio  = a = 20 

(c)  Initial  temperature  of  the  mixture  = 1200°  R 
[This  corresponds  to  a compression  ratio  of  8 
for  adiabatic  transformation] . 

(d)  Changes  of  pressure  and  specific  heat  ratio  across 
the  flame  front  are  neglected. 

(e)  The  motion  of  the  power  piston  during  the  combustion 
processes  is  negligible. 

Construction  of  the  wave  diagram  is  based  on  standard 
methods  given  in  references  (2)  and  (3).  The  movement  of  the 
interface  with  respect  to  time  as  combustion  proceeds  to  com- 
pletion is  shown  by  the  dotted  lines  on  the  wave  diagram. 

After  a certain  time,  when  the  flame  front  reaches  the  inter- 
face, combustion  stops. 


ANALYSIS  OF  THE  WAVE  DIAGRAM 


It  is  interesting  to  note  that  the  strength  of  the 
initial  shock  wave  characterized  by  the  shock  mach  number  is 
comparatively  low  (Mg  = 1.247)  for  the  assumed  typical 
initial  conditions.  The  progress  of  the  computations  is 
recorded  in  Table  1 which  give  the  Riemann  parameters  and  the 
pressure  ratios  computed  for  the  various  regions  labelled  in 
the  wave  diagram  (Fig.  4). 

It  is  clearly  seen  from  Figure  4 that  once  the  interface 
reaches  the  balancing  chamber,  the  reflected  waves  from 
previous  interactions  moving  in  the  opposite  direction  towards 
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the  main  chamber  slows  down  the  interface  and  after  some  time 
the  interface  itself  starts  moving  towards  the  main  chamber. 
Because  the  reflected  shock  waves  are  all  of  comparatively 
small  mach  numbers  in  the  neighborhood  of  1.1  the  effect  of 
the  interaction  of  these  waves  on  the  flame  front  itself  is 
negligible.  However  towards  the  end  of  the  combustion  process 
the  flame  front  is  appreciably  slowed  down  and  finally  when 
the  flame  front  reaches  the  interface,  the  combustion  stops. 

For  this  specific  case  the  combustion  time  is  found  to  be 
of  the  order  of  0.15  millisecond  which  is  one  order  of  magni- 
tude lower  than  that  observed  from  the  actual  engine  pressure  - 
volume  indicator  diagram.  This  discrepancy  is  most  likely  due 
to  assumption  (e)  where  the  motion  of  the  power  piston  during 
combustion  was  neglected  whereas  in  the  actual  more  general 
situation  this  motion  cannot  be  neglected. 

It  should  be  noted  that  in  the  absence  of  any  information 
pertaining  to  the  non-steady  flow  process  in  the  NAHBE  engine, 
construction  of  the  wave  diagram  with  the  simplified 
assumptions  mentioned  above  as  a first  step  give  valuable  in- 
sight into  the  flow  dynamics.  Further,  the  information  thus 
obtained  will  in  turn  help  refine  the  analytical  model  for 
subsequent  analysis. 
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DISCUSSION  OF  HEAT  RELEASE  RATE 

The  rate  of  heat  release  of  a burning  gas  mixture  as  a 
function  of  the  flow  conditions  is  very  much  dependent  on  the 
particular  combustion  process  and  flow  geometry.  Burning  pro- 
ceeds throughout  the  mixture  in  a complicated  manner.  It  thus 
becomes  necessary  to  select  some  heat  release  model  on  which 
the  wave  diagram  procedures  can  be  based.  As  *n  example,  one 
may  assume  for  instance  a linear  heat  release,  namely, 

Q - Qt  (~t  — ) ----  (2) 

c 

where  t.  is  the  initial  (or  ignition)  time  and  t the  total 
duration  of  the  combustion  process.* 

The  pressure  equation  thus  obtained  based  on  the  linear 
heat  release  model  is  repeated  here  for  convenience.  Thus  we 
have 

cic2  * ^ I**  ^ * vr  (a_b)  tn  ivrwt  * /v?i 


3/2 

(vt) 


(3) 


* The  duration  of  combustion  in  the  following  analytical 
development  for  numerical  calculations  is  taken  to  be 
70  degrees  crank  angle. 
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The  constant  c.^  in  equation  (3)  is 
pressure  in  the  combustion  chamber  at  the 
the  total  heat  of  reaction,  and  the  total 
The  motion  of  the  piston  is  given  in 
Volume,  Vt , by  the  expression 


determined  by  the 
time  of  ignition, 
time  of  combustion, 
terms  of  the  total 


V = A - B cos  (wt)  (4) 


1.1 

ll1 


where  A,  B,  and  w are  constants  to  be  evaluated  for  the  given 
engine  geometry. 

Now  let  0 - wt  (5)  be  the  crank  angle  from  the  top  dead 
center.  Thus  |a  (0)  can  be  calculated  from  equations  (3)  and  (4). 

For  the  particular  engine  geometry  under  investigation  in 
the  present  study,  the  value  of  the  constants  A and  B are 
given  by 


A = 

27.22 

B = 

22.05 

- 20° 

< 0 < 20 

A = 

25.3 

B = 

20.1  , 

O 

o 

C\j 

A 

a> 

v 

i A 

to 

o 

At  high  loads  the  (0)  history  during  combustion  calcu- 
lated from  equation  (3)  for  y • 1.5  give  higher  pressures  com- 
pared to  the  experimentally  observed  data  particularly  later 
in  the  expansion  stroke.  This  is  shown  in  Figure  6 which  also 
shows  the  effect  of  specific  heat  ratio,  y,  on  the  calculations. 

For  y - 1.35  the  calculated  pressure  differ  from  the  experi- 
mental values  during  the  initial  phase  following  ignition  only 
slightly.  However,  considerable  deviation,  though  not  as  pronounced  as 
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for  y = 1.5,  is  still  observed  later  in  the  combustion  cycle. 

In  order  to  improve  on  this  heat  release  model  from  the 
point  of  view  of  matching  theoretical  and  NAHBE  experimental 
results  an  attempt  was  made  to  change  the  nature  of  heat  re- 
lease to  a form  more  compatible  with  that  observed  from  actual 
spark  ignition  engine  tests.  Data  from  several  engine  tests 
and  reported  in  reference  (4)  have  shown  that  the  actual  heat 
release  diagram  follow  a "S"  shaped  curve  as  shown  in  Figure  5. 

Therefore  a family  of  different  heat  release  curves  having 
the  general  form 

- a f (0)  ----  (6) 

where  a = constant 

were  examined  and  subsequently  used  in  the  derivation  of  the 
pressure  equation.  The  pressure  levels  computed  from  this 
equation  were  compared  with  the  experimental  NAHBE  test  results. 
On  the  basis  of  this  comparison  the  following  equation  for  the 
heat  release  model  was  used  for  subsequent  analysis. 


Where  Vt  = total  Volume. 

Equation  (7)  gives  a simple  expression  for  pressure 
(Equation  8)  in  terms  of  two  constants  a and  Cj^  the  values 
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of  which  can  be  evaluated  from  two  known  p (0)  values. 
(y-1)  a 0 ♦ cic2 


This  is  most  easily  accomplished  for  computational  purposes 
from  the  engine  indicator  diagram  using  two  "arbitrary"  points 
namely  ignition  point  and,  say,  30°  ATDC. 

It  has  been  found  necessary  from  experience  that  these  two 
"arbitrary"  points  should  be  located  sufficiently  well  apart 
for  optimum  results. 

Figures  (6) , (7) , and  (8)  compare  the  results  based  on  the 
two  heat  release  rates  with  NAHBE  experimental  results  at  1500  rpm. 
The  following  conclusions  can  be  drawn  from  this  analysis: 

(1)  At  high  loads  the  he3t  release  rate  equation  (7)  with 
Y = 1.35  give  good  agreement  with  the  experimental 
data  during  the  initial  phase  of  combustion  over  a 
period  of  30°  crank  angles  from  the  ignition  time. 

However  the  computed  pressures  are  found  to  be  20  to 
80%  higher  than  the  experimental  values  later  in  the 
combustion  cycle. 

(2)  At  very  light  loads  the  two  heat  release  rates  give 
comparatively  good  correlation  with  NAHBE  experimental 
data  over  the  entire  combustion  phase. 

It  is  of  considerable  interest  to  compare  the  above  results 
with  experimental  data  from  a conventional  otto  cycle  engine 
for  identical  test  conditions.  Since  accurate  pressure  traces 
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from  the  latter  are  not  available  at  this  time  it  is  recom- 
mended that  such  a comparison  be  made  as  otto  cycle  engine 
test  data  becomes  available.  This  could  conceivably  shed 
some  light  as  to  the  nature  of  NAHBE  combustion  process. 

VISUAL  STUDY  OF  NAHBE  COMBUSTION  AT  VERY  LIGHT  LOADS 

Direct  visual  observation  through  high  speed  photo- 
graphs of  NAHBE  combustion  process  with  ethyl  alcohol  as  fuel 
at  very  light  loads  conducted  on  an  experimental  square  engine 
at  1000  rpm.  have  revealed  some  interesting  facts. 

1.  A reaction  zone  of  relatively  large  size  and  of  pale 
blue  color  is  observed  immediately  after  ignition 
followed  by  multiple  hot  spots  (glows). 

2.  In  certain  frames  the  combustion  intensity  increases 
towards  the  middle  of  the  expansion  stroke. 

3.  Glowing  of  particles  towards  the  end  of  the  expansion 
stroke  and  a subsequent  sweeping  action  even  while  the 
piston  moves  down  is  generally  observed  in  the  main 
chamber.  The  same  phenomena  also  occurs  at  times  in 
the  balancing  chamber. 

4.  Appearance  of  a blue  region  even  before  ignition  in 
certain  frames  suggests  a long  ignition  delay. 

5.  The  combustion  duration  appears  to  be  considerably 
larger  than  that  originally  anticipated. 
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l 4 


l 

While  these  visual  observations  are  of  considerable 
interest,  it  should  be  pointed  out  that  the  photographs  were 
taken  at  a relatively  slow  speed  of  200-250  frames  per  second 
and  at  very  light  loads  because  of  practical  considerations. 

In  order  to  gain  further  understanding  of  the  NAHBE  combustion 
process  a more  realistic  approach  would  be  that  of  increasing 
the  frame  speed  (which  helps  observe  the  particle  path  lines) 
and  conducting  some  tests  under  moderate  load  conditions.  Thus 
the  combustion  phenomena  at  different  loads  can  be  studied. 
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FIG.  1.  ANALYTICAL  MODEL  FOR  VOLUME  MODE  ANALYSIS 


MODE  U NOMINAL  SPEED  = 150C  RPM 


DYNAMOMETER  LOAD  - 9.562  FTLBS 


10°  □ Linear  Curve;  rj  = 100%  = 65 

7 =1.35 


© NAHBE  Experimental 


■ Linear  Curve  ; rj  = 100% 


TOTAL  VOLUME,  CUBIC  INCHES 


FIGURE  6 


Nominal  Speed  1500  RPM 
MODE  R 


DYNAMOMETER  LOAD  - 7.562  FTLBS 


□ Linear  Curve  ; r}=  100%  , 9f 
7=  1.35 

© NAHBE  Experimental 


TOTAL  VOLUME,  CUBIC  INCHES 


FIGURE -7 


MODES 


DYNAMOMETER  LOAD  - 6.26  FTLBS 


Nominal  Spsad  1500  RPM 


O Linear  Curve 
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TOTAL  VOLUME,  CUBIC  INCHES 


FIGURE  -8 


Table 

1.  Data 

Sheet  for 

Nonsteady 

Flow  Analysis. 

Region 

P 

Q 

U 

A 

S 

VK 

1 

5.0 

5.0 

0 

1 

0 

1 

2 

5.75 

5.0 

.37 

1.076 

.009 

1.65 

3 

5.75 

5.0 

.37 

1.076 

.009 

1.65 

4 

5.75 

5.3705 

.19 

1.112 

.010 

2.08725 

5 

5.95 

5.0 

.475 

1.095 

.017 

1.86 

6 

5.95 

5.1765 

. 38675 

1.11265 

.017 

2.083 

7 

6.  02 

5.0 

. 51 

1.102 

.017 

1.94 

8 

6.02 

5.135 

.4425 

1.1155 

.017 

2.1146 

9 

6.07 

5.0 

. 535 

1. 107 

.024 

2.0 

10 

6.09 

6.09 

0.0 

1.21825 

.0435 

3.81 

11 

6.0295 

6.0905 

-.0305 

1.212 

.0435 

12 

6.02 

6.0802 

- .0301 

1.210 

.029 

3.695 

13 

-.2394 

2.0425 

3.096 

1.947 

14 

5. 6395 

5.3705 

.1345 

1.  101 

.010 

1.947 

15 

10.4525 

10.4525 

0 

2.0905 

3.096 

2.2877 

16 

10. 452 

10.598 

- .073 

2.105 

3.096 

2.402 

17 

5.97 

5.3704 

.29965 

1.134 

.011 

2.402 

18 

5.97 

5.4924 

. 2388 

1.14624 

.011 

2.582 

19 

6.215 

5.00 

. 6075 

1.1215 

.017 

2.1985 

20 

10.60 

10.60 

0 

2.12 

3.096 

2.5221 

21 

10.453 

10.746 

-.1467 

2.12 

3.096 

2.5221 

22 

5.9316 

5.4924 

.2196 

1.1424 

.011 

2.5221 

23 

6.215 

5.1765 

.51925 

1. 13915 

.017 

2.46 

24 

6.  215 

6.069 

.073 

1.228 

.028 

4.10 

n * 


25 

10.453 

10. 297 

.078 

2.075 

3.096 

2.169 

26 

10. 546 

10.192 

.177 

2.0738 

3.096 

2.160 

27 

6.  1869 

5.0 

. 5934 

1.1187 

.017 

2.16 

28 

10.5483 

10.4017 

.0733 

2.095 

3.096 

2.31 

29 

6. 1137 

5.1763 

. 4687 

1.129 

.017 

2.31 

30 

10.5513 

11.5487 

-.4987 

2.21 

3.10 

3.3495 

31 

5.8612 

6.069 

-.1038 

1.193 

.028 

3. 3495 

32 

6.215 

6.06 

.0775 

1.2275 

.028 

4.088 

33 

6.1864 

6.08 

.0532 

1.2266 

.029 

4.0774 

34 

6.185 

6.185 

0 

1.237 

.029 

4.322 

35 

10.60 

11.2095 

-.30475 

2.181 

3.097 

3.070 

36 

11.21 

11.21 

0 

2.242 

3.098 

3.7147 

The  relatively  large  number  of  wave  passes  per  cycle 
appears  to  justify  a "volume  mode"  approach  (whereby  the 
pressure  at  each  instant  is  assumed  to  be  uniform  throughout) 
as  a first  step  in  the  analysis.  What  is  sought  here  is 
the  effect  of  pressure  exchange  on  NAHBE  performance.  The 
error  that  may  result  from  the  quasi-static  treatment  of 
pressure  exchange  that  is  implied  in  the  volume-mode  approx- 
imation is  probably  insignificant,  because  shock  waves  are 
unlikely  to  form  within  the  limited  space  of  the  NAHBE  chambers 
anyway.  More  serious  error,  on  the  other  hand,  will  surely 
result  from  three  assumptions  that  are  dictated  here  both  by 
the  need  to  make  the  model  analytically  manageable  and  by 
the  lack  of  information  about  the  details  of  the  actual  process 
These  assumptions  are: 

(1)  that  no  mass  or  energy  is  transported  across  the 
the  interface  between  the  combustion  gas  and  the 
"balancing  gas"; 

(2)  that  flow  losses  can  be  neglected,  and 

(3)  that  all  the  transformations  of  the  balancing  gas  are 
isentropic. 

Fig.  1 describes  the  analytical  model  used  in  the  present 
analysis.  Region  1 is  occupied  by  the  combustion  gas,  region 
2 by  the  balancing  gas.  The  interface  between  the  two  is  simu- 
lated by  a massless,  frictionless,  impermeable,  and  adiabatic 


piston.  This  piston  is  free  to  float,  thus  allowing  the  pres- 
sure to  be  equalized  throughout  regions  1 and  2 at  each 
instant. 


Fig.  1 


The  following  symbols  will  be  used: 

c specific  heats  at  const,  pressure,  const,  volume, 

respectively 

mass  of  gas 

pressure 

heat  released  in  the  combustion  process 

specific,  total  entropy 

time 

static  temperature 


specific,  total  volume 


V°v 


Subscripts  1 and  2 will  refer  to  regions  1 and  2,  respec- 
tively; subscripts  i and  f will  denote  inital  and  final  con- 
ditions, respectively,  and  subscript  t will  denote  'tota it 
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The  power  piston  will  generally  be  moving  during  the 
combustion  process.  The  effect  of  this  motion  will  be  fully 
accounted  for  in  the  present  analysis.  However,  to  acquire 
an  immediate  "feel"  for  the  effect  of  pressure  exchange,  con- 
sider first  the  case  in  which  the  power  piston  is  either 
stationary  or  moving  so  slowly  that  changes  of  the  total 
volume  vt=vi+v2  dur*n9  the  combustion  process  can  be  neglected. 
Under  these  conditions, 

dQ  * *C„  0*i  ) 

« (pdV,  ♦ Vf  <kjo  + 

=r  J~  \4 

r-i 

• • 

p.  - p;  = (V-  I ) Qt 

independently  of  the  presence  or  absence  of  a balancing  chamber. 

Pressure  exchange  has,  however,  a significant  effect  on 
the  overall  entropy  production.  Indeed,  by  virtue  of  assumption 
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and 


AS  = m,  A*, 


This  quantity  is  always  negative  because,  9^  being  greater 
than  1.0,  + all  positive  terms.  Therefore,  the  proposed 

utililization  of  pressure  exchange  has  the  effect  of  increasing 
the  available  work  for  any  given  energy  input  and  total  volume. 
Note,  however,  that  this  effect  is  independent  of  the  initial 
temperature  in  the  balancing  chamber. 

Consider  now  the  more  general  situation  — that  in  which 
the  motion  of  the  power  piston  during  the  combustion  process 
is  not  negligible.  Let  this  motion  be  a harmonic  one,  such 
that 

\f^.  5 =*  A - j?>  'Cos.  cut" 

where  A, B , and  u are  constants,  and  let 

where  t^  is  the  initial  (or  ignition) time  and  tc  the  total 
duration  of  the  combustion  process.  Then, 

\4  = A - <u>*  ^ 


4 


where 


i - *“  Q ♦ * ) 


5 


and 


U ! 

toXc.  A 


From  (2) , 


« 

u (a  ~ f~>  £*3  ^ ^ = c,  ( 


a constant) 


Substituting  into  Eq.  (3) , we  obtain 


“***  Cl  0 la  - Wfc 

CV-.J  <Sft  * 

This  equation  may  be  integrated  in  closed  form  for  f = */£. 
and  through  the  use  of  the  relation  cos^=  I which  provides 

a good  approximation  over  the  range  of  ^ within  which  combustion 
normally  takes  place  in  reciprocating  engines. 


Then, 


C,  1-  / <s, 


where  Cj  is  another  constant.  Finally,  remembering  that  p = uz  and 


that  |- -£-=coswt. 


CiCt  f gjL.Ltfa-B'+rf ■ * 


A — 3 c-o*  u>r 


The  product  C^C2  now  forms  a single  constant,  the  value  of 
which  is  determined,  through  Eq.  (6) , by  the  position  and  velocity 
of  the  piston  at  ignition  time  (t=ti) , the  total  heat  of  reaction, 
the  engine  rpm,  the  total  time  of  combustion,  the  volumetric 
compression  ratio,  and  the  pressure  in  the  combustion  chamber  at 
the  time  of  ignition. 

The  average  i”  in  combustion  is  actually  around  4/3.  The 
use  of  f = 1.5,  as  has  been  done  above  for  the  sake  of  analytical 
simplicity,  introduces,  therefore,  a quantitative  error  in  the 
results.  It  does  not,  however,  affect  the  validity  of  the  basic 
conclusion  that  may  be  drawn  form  Eq.  (6),  namely,  that  the  p(t) 
history  during  combustion  is  independent  of  the  presence  or  ab- 
sence of  a balancing  chamber,  just  as  in  the  case,  previously 
treated,  of  the  stationary  piston.  This  result  could,  of  course, 
have  been  predicted  on  the  basis  of  the  observation  that  the 
volume  fraction  <f>  does  not  appear  in  Eq.  (1)  . 

The  entropy  rise  is  again  obtained  as 

AS  * m,  (<v  u h + e*  ) 
v A f ' 


In  this  case,  pi  and  ti  being  given,  pf  is  calculated  from  Eq.  (6) 
for  t = t.  + t . Then, 

- * (ff 

hence , . . 


Vfi  = vtt  - 14  (Jg-) 

-tofo-YCtf] 
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Therefore, 


Now, 


antilogarithm  in  Eq.  (7)  is  always  greater  than  1.0,  and  it 
follows,  as  in  the  case  of  the  stationary  piston,  that  the  proposed 
utilization  of  pressure  exchange  has  the  effect  of  increasing 
the  available  work  for  any  given  energy  input. 


* It  would  be  1.0  only  if  the  entire  process,  in  both  chambers  1 and 
2,  were  isentropic. 
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APPENDIX-2 


Derivation  of  Equations  and  Sample  Calculations:  - 


Nomenclature:— 


Speed  of  Sound  (dimensional) 

Speed  of  Sound  (Nondimensional) 

Heat  Input  per  unit  weight  of  Fuel 

Mechanical  Equivalent  of  Heat 

Constants 

Mach  Number 

Pressure 

Riemann  Variables 

Heat  released  in  the  Combustion  process 
Entropy  (Nondimensional) 

Temperature 
Velocity  (dimensional) 

Velocity  (Nondimensional) 

Specific  Heat  Ratio 
Cross-Sectional  Area;  Air-Fuel  Ratio 


Subscripts 


Burned  Condition 

Combustion 

Flame 

Initial 

Shock 

Unburned  Condition 


D 

0 


i 


I 


i 


Definition  of  Riemann  Variables:  - 

The  Riemann  Variable*  are  defined  by  the  following  equation*. 

P = y~T~  A + U - la) 


Q = j ■ A - U -(b) 


For  7 - 1.4  we  then  have 
P = 5 A+U  -(c) 

0 = 6 A-U  -(d) 

From  (c)  and  (d)  we  get 


U = 


2 


and  A = 


10 


Nonsteady  Flow  Analysts  of  NAHB  Combustion  Process: — 

Derivation  of  Equations  and  Sample  Calculations: 

Consider  the  flame  propagation  in  a constant  area  duct  as  shown  in  figure  (1)  below. 


Ub 


Burned  Gas 


T 


Vf 


Uu 


Unburned  Gas 


Flame  propagating  into  unburned  gas 
moving  at  Velocity  Uu 


Constant  area  duct 

Fig.  1.  Flame  propagation  in  a constant  area  duct. 


From  the  conservation  of  mass  and  energy  and  from  the  equation  of  state  the  following 
equation  can  be  derived. 


Ub  = Uu  - V\ 


(1  + 


0 

0 

0 

fl 


II 

0 


2 


t 


u 


a 


4# 


D 


D 

0 


Equation  VI.  h.  8 of  reference  (3) 
where 

K,  = (7U“  D 9 J Q/a0J 


For  the  case  of  constant  specific  heats  on  both  sides  of  the  flame  front  we  have 

yu  = ?b 

and  equation  ( 1 ) reduces  to 

K, 


Ub  = Uu  - Vf 


A 2 
Mu 


(2) 


Furthermore,  assuming  that  the  burning  velocity,  vf  is  proportional  to  the 
temperature  of  the  unburned  gas. 


vf  = K Tu  = K 


vi  - 

' 1 8o 


au2 


T'u  Rg 
a0  Au2 

tjtt 


(3) 

(4) 


Substituting  equation  (4)  into  equation  (2)  we  get 


Ub  = Ug  - = 


KHtjc 


(5) 


K Q = KHt?c 
Cp  aQ  Cp  a0  (a  + i ) 


(6) 


cp  a°  Cp  a0  (a  + 1 ) 

However,  since  the  burned  gas  is  adjacent  to  a closed  end  of  the  duct,  we  have 
Ub  = 0 

and  from  equation  (5),  Uu 
where 

Q = Heat  Input  per  unit  weight  of  fuel. 

_ Ht?c 

(«+  1) 

Let  the  initial  conditions  of  the  gas  in  the  duct  be  such  that: 

H = Heating  value  of  fuel  = 19,000 
a = Air-Fuel  ratio  = 20 

T = 540°  R 


Btu 


'°0  ss 


r?c  = combustion  efficiency  = 90% 

To  = Initial  temperature  of  air  = 1200°R 

ft 


a0 
Then  v 


Initial  speed  of  sound 
ft 


1700 


sec 


f 


222.22  @ 1200  R 


< - - 0 ,85 


From  equation  (6),  Uu  = 0.370 

The  gas  ahead  of  the  flame  front  is  accelerated  from  rest  to  this  velocity. 
The  strength  of  the  Initial  Shock  is  given  by 

I A u I Uu  - U0 


A A0 

Then  we  have  from  shock  tables  1 a of  reference  (3) 


= 0.370 


= Au  = 1.076  ; f-  = 1.65 

Ms  = 1 .247  •*=  Shock  Mach  Number 

AS  = Su  - S0  = 0.009 

v Ys  222.22 

» a0  1700 


= 0.13 


Wf  =■  Velocity  of  flame  front  relative  to  the  duct 
= Uu  + K,  = 0.37  + 0.13  = 0.50 


Also,  K,  = ( 7 ~ ap29  J ° = 2.8234 
Ab  = Au  (1  + ■—-)  x'2  = 1.995 

MU 

Sb  = su  + Cn 

1.995 

= 0.009  + 5 8n  ( ) = 3.096 


Thus  for  the  above  interaction,  the  numerical  results  obtained  are  shown  in 
Figure  (2)  Below. 


1 


0 

0 

0 

0 

n 


1 

• - 

u 


Fig.  2.  Flame  Front  advancing  into  a combustible  mixture. 
(Ignition  at  a closed  end  of  the  duct) 


B 

n 


n 

0 

1 
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Interaction  of  the  Initial  Shock  (M,  = 1.247)  at  the  first  change  of  cross-section: - 

Let  subscripts  L and  R denote  conditions  immediately  to  the  left  and  to  the  right, 
respectively,  of  a discontinuous  are  change.  Assuming  that  the  flow  is  subsonic  on 
both  sides,  the  step  section  is  reached  by  an  p characteristics  from  the  left  and  by 
a Q characteristics  from  the  right.  Therefore,  PL  and  QB  are  known. 


Fig.  3.  Shock  at  a discontinuous  change  of  cross-section. 


Since  the  shock  is  weak,  P4  = P3  = 5.75 
and  Qs  = Q,  = 5.0 

Then  for  regions  (4)  and  (5) 


5.75 

5.0 


1.15 


0(R 

and  = 0.44 

«L 


With  the  above  known  values  go  to  figure 

6e  - 4 (FOA).  We  then  have 
Ql  Q4 

o—  = -5-  = 0.934  - Q4  = 5.3705 

' L “4 

and  -77-  = 77-  = 1.19  - P5  = 5.95 

Mr  Mj 


A4  = 


P4  + Q4 

10 


1.112 


u4 


P4-Q4 

2 


0.19 


.i 


5 


i 


(1)  ( A Q)n  _ 0.3705 

A3  1.076 


= 0.344.  Hence  M,n  = 1.1085 


AS  - 0.001  -g-  = 1.265 


p3  = 1.65  p0  p4  = 1-265  (1.65  p0) 
= 2.08725  p0 


Thus  the  reflected  shock  is  of  Mach  number  1.1085. 

(2) 


( A P ),'  _ 0.95 


= 0.95  . Hence  M,  = 1.318 


^ 10  U,9W  • I IOI Iflj 

AS  - 0.017  y = 1.86 

ps  = 1.86  p0 

Thus  the  transmitted  shock  is  stronger  than  the  incident  shock. 


The  above  computations  are  summarized  below  in  tabular  form. 


p = 

2 

7-1 

- A - u 

A = 

P + Q 
10 

Q = 

CM  1 

- A - u 

U = 

P-Q 

2 

Region 

P 

Q 

U 

A 

S 

P'/Po 

1 

5 

5 

0 

1 

0 

1.00 

2 

5.75 

5 

0.37 

1.076 

0.009 

1.65 

3 

5.75 

5 

0.37 

1.076 

0.009 

1.65 

4 

5.75 

5.3705 

0.19 

1.112 

0.01 

2.08725 

5 

5.95 

5.0 

0.475 

1.095 

0.017 

1.86 

Following  the  method  described  above,  the  nature  of  the  subsequent  interactions  of 
the  transmitted  shock  (M,  * 1.318)  and  the  resulting  transients  at  downstream  discon- 
tinuous change  of  cross-sections  can  be  determined.  The  details  of  the  computations  are 
given  below. 


6 


0 


Fig.  4.  Shock  at  discontinuous  change  of  cross-section. 

Case  1. 

Case  2. 

S--M 

It  = 0.85 

Since  the  incident  shock  is  weak, 

Since  the  incident  shock  is  weak, 

P6  = Ps  = 5.95 

P,  = P7  = 6.02 

and  Q7  = Q,  = 5.0 

and  Q9  = Q,  = 5.0 

Then  for  regions  (6)  and  (7) 

Then  for  regions  (8)  and  (9) 

ofe"  = (57  = 1,19  and  = 08 

= ■?*-  = 1.204 

Qr  Q» 

With  the  above  known  values  to  figure 

6 e - 4(FOA). 

With  the  above  known  values  from  figure 
6e-4 

Then  from  figure 

-2k  » -9l.  = 0853 

- y5-  = 0.870 

pL  p6 

- Q6  = 5.1765 

PL  P,  0853 

■*  Q,  = 5.135 

Case  1,  continued. 


Case  2,  continued. 


Shock  @ closed  end:— 


Strength  of  reflected  shock 


= 0.535 
1.107 

=»  Ms  = 1.331 

_Al 

A 


0.483  and  U10  = 0 

l 

— = 1.905 
P 

1.1005  - Al0  = 1.21825. 


P,  = 2 Po 


AS  = 0.0195  and  Q, 0 = P| o = 5 A, 0 = 6.09 
p10  = 1.905  (p,)  = 3.81  p0 


Reflected  shock  at  change  of  cross-section. 


Fig.  5.  Shock  at  change  of  cross-section. 

Q,o  =Q.i  = 6.09)  PL  _ Pia  _ 6.02  _ _ qooc 

Ps  = P, , = 6.02  | Qr  Q, , 6.09  0 9885 

P P 

From  figure  6 e - 7 ; = —■  = 0.99  P, , = 6.029 

Qr  Qi  i 

Qi  Qi  2 

-pj  = * 1.01  Q,2  = 6.0802 

9 


J 


= 0.8473  ■*  Ms  = 1.282 


AS  = .012 
-Sr  = 17475 


Pii  -Pio  _ 6.029  - 6.09 
A10  1.21825 


EXPN  WAVE 


Ui,  = (P,,  + Q,  ,)/2  = -0.0305 


Hence  the  flow  in  region  (11)  has  a negative  velocity. 


Ai  i 


(P.i  - Qm) 
10 


= 1.212 


Similarly  for  region  12  , we  have 

P,2  = 6.02  and  Q,j  = 6.0802 
Hence  U12  = -0.0301  and  A,,  = 1.210 
S,j  = S8  + 0.012  = 0.033 
p,2  = 1.7475  (p„)  = 3.695  p0 


Thus  all  flow  properties  in  regions  (11)  and  (12)  are  completely  defined. 

Note  that  the  transmitted  shock  because  of  passage  into  a larger  cross-sectional  area. 
Thus  this  interaction  results  in  a weak  shock  and  an  expansion  wave. 


Merging  of  shocks:— 

With  reference  to  figure  6g  - 4 of  (FOA)  we  have  the  following  situation. 


II 


II 

0 

0 

Fl 


0 


r 


11 


li 


n 


11 


u 


i. 

0 

0 


For  the  case  of  weak  shocks,  we  have 

(AAlm  = (A  A),  + (AA)„  = 0.108 
Since  (A  A),  = A8  - A7  = 0.0135 
and  (AA)n  = Ai2  - Ag  = 0.0945 


Then  we  get 


(AA),,,  = An  “ A7  = 0.108 


J — i i 0-108  _ f nqo 

A,  “ 1 + 1.102  1 098 

p' 

Now  froM  shock  tables  we  have  -jj-  = 1 .87  and 

Ms  = 1 .322  propagating  into  region  7 . 

A,,  = 1.098  A,  = 1.210. 

It  can  be  seen  that  the  merger  produces  a stronger  shock. 

Interaction  of  Flame  Front  with  the  First  Reflected  Shock: 

Figure  (7)  below  shows  a typical  interaction  involving  a flame  front  and  a shock 
wave.  The  matching  condition  across  the  flame  front  are  the  pressure  and  flow  velocity. 


Guess  A,  3 = 2.1 


= 1.0526 


Now  from  shock  tables  we  have 


M.  = 1.1685 


f = 1.425 


0.2615  AS  = .0035  etc. 


U13  = — 0.5217 


Note:- 


For  a Q shock 


I A u I U ahead  — U behind 

A A a head 


Hence  0.2615  = - rU>3 

Ab 

and  since  U5  = 0 (Velocity  adjacent  to  the  closed  wall) 
U13  = - 0.2615  (1.995)  = - 0.5217  «= 


Also  = [1.425(1.650)  = 2.35125] 

Hence  p13  = 2.35125  p,  = p,4  [Matching  Condition] 


’ten  -jjj1-4  = 1.1264  =■*■  Ms  = 1.0505  (From  shock  Table) 
A S = o 


= 0.0825 


-^  = 1.0165 


U,  4 - U4 

^ = 0.0825  This  assumes  that  the  resulting  transient  is  a shock 

U,«  = U«  + 0.0825  (1.112) 

U,  ^ = 0.2817 

Hence  the  difference  in  flow  velocity  across  the  flame  front  after  the  interaction  is  given  by 
AU  = U14  - U,j  = 0.09826- (-0.5217) 

= 0.80344*  0.37 

Hence  solution  is  Not  correct 
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2nd  Trial  (Identical  to  the  procedure  describe  above) 


Guess  A,  3 = 2.0425 


Then  -vL2=  1.0238 
Ab 


From  shock  table,  we  have  the  following:— 


Ms  = 1.074  AS  = 0 -g- 


iMLL  . o 120  = 


Ub  = o 


Ub  - U13  = 0.2394  -*  U,3  = - 0.2394 
p,  3 = 1.947  p,  = P14  (Matching  condition] 


^ = 0.9328  = (41-2-) 
m4  m4 


For  Isentropic  Waves  since  p,  4 < p4 


Thus  the  resulting  transient  is  an  expansion  Wave.  And  for  a P - expansion  Wave 
Q,4  = Q4  = 5.3705 


= 0.990 
A4 


a14  = 1.101 

U14  = 5 A|4  — Q|4 
U14  = 0.1345 


AU  = U14  - U,j  = 0.1345  + 0.2394  = 0.3739 

Velocity  t L Velocity  behind 
ahead  of  the  flame-front 

flame-front 

Hence  A u]  = 0.37  as  0.3739  * A U]  *■  co^ditiwi 


before 

interaction 


Thus  tho  trial  solution  is  Correct 


After 

Interaction 
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Shock  at  closed  End:— 


Hence  A,  , = 2.0905  p,  s - 2.2877  p, 

A S = 0 ■*  S,j=S13  = 3.096 

The  Reimann  parameters  are  determined  as  follows:— 

P,  s = 5 A,j  + U,s  = 10.4525 
Qis  = 5 A,,  ' UIS  = 10.4525 
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Shock  at  change  of  cross-section:— 


L R 


Fig.  10.  Shock  at  Area-change 


Since  the  Incident  shock  is  very  weak  we  have 

Pi  7 = Pi  8 = 5.97  and 
Qi  9 = Qs  = 5.0 

Hence  - u94- 

p 

With  the  known  value  of  z: — go  to  figure  6 e - 4 and  obtain 

Ur 

IT-  = = 0.92  =»  Q,  8 = 5.4924 

rL  rl  8 

and  -fe = fe" = 1243  ■*  p«’ = 6215 

U,»  = 0.6075  and  Al9  = 1.1215 
U | § = 0.2388  and  Al8  = 1.14624 

Now  for  the  transmitted  shock 

LA_p  = fis ,~_ps  = 0242  - Ms  = 1.0746  and  = 1.182 


And  for  the  reflected  shock 


(A  Q)  Qis  Q17  p' 

— r-  = — : = 0.1076  - Ms  = 1.033  and  y = 1.075. 
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Crossing  of  a contact  surface  by  a shock  Wave:— 


i 


I 


TT 
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This  type  of  interaction  appears  in  a wave  diagram  as  shown  in  figure  (11)  for  a Q - shock 
crossing  an  interface 


/ 

/ 


Fig.  11.  Shock  Wave  crossing  an  Interface. 


It  is  to  be  noted  that  for  a weak  Q - shock  I crossing  an  interface  with  7R  = 7,,  in 
general  we  may  expect  a transmitted  wave  I'  and  a reflected  wave  II  to  emerge  from 
the  crossing  point.  However  in  the  present  case,  the  interface  which  merely  separates 
the  combustible  mixture  of  gases  from  air,  is  unaffected  by  the  passage  of  the  shock 
because  of  the  absence  of  any  entropy  gradients  across  the  interface.  However  the 
interface  is  slowed  down  by  the  incident  shock  as  can  be  seen  by  comparing  the  flow 
velocities  across  the  interface  in  regions  (5)  and  (6)  respectively. 


Thus,  we  have  Us 
and  U6 


0.475 

0.38675. 
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Flame  Front  at  First  Change  of  Cron-section: 


Fig.  12.  Flame  Front  6 change  of  cron-section 


The  matching  condition  acron  the  flame  front  are  the  pressure  and  A U.  Assuming 
that  the  resulting  transient  is  an  Isentropic  expansion  wave  acron  regions  (21)  and 
(25)  we  have 

PJS  = P, , = 10.453 


Now  Guess  A 


) for  an  Isentropic  Wave 


0.0376 


0.0376  (From  the  table  of  Mach  Number  Functions) 


Now  from  the  conservation  of  mass  for  regions  (25)  and  (26)  we  get 
Dj  $ fljj  = D2  $ o j 6 

D = D Uls  = d rtt|-  ^ = °2i 

°2‘  °2S  026  °2S  (oR  ’ 0.44 


Dj6  = 0.0852  -*  M26  = 0.0855  (From  the  table  of  Mach  Number  Functions) 
From  the  conservation  of  Energy  for  regions  (25)  and  (26)  we  have 

A2J2  + Uj,2  = A262  + ^ U262 

or  5AJS2  + U252  = 5 A262  [1  + 0.2M26j] 

Hence  21.5340  = 5.00731  A262 
A26  = 2.07377 
and  U16  = 0.177 

= 0.860096  - pJ6  = 2.169  p, 


Now  for  the  Isentropic  flow  regions  (25)  and  (26)  we  have 

sir  - <-b’=  - 099585 


Hence 


p26  = 2.169  (0.99585)  p,  = 2.16  p, 
p27  = p26  - 2.16  p,  (Matching condition) 


r — = 0.9825  Expansion  Wave 

9 


Hence 


-r12-  = (0.9825)  V7  = 0.99748 

At  9 

A,  7 = 1.11867 

Qj  7 = Qi  9 - 5.00  across  the  expansion  Wave 
U27  = 5A27  - Qj,  = 0.59337 
A U = U2  7 - U26  = 0.5934  - 0.177  = 0.416 


AU]  - 0.37  a 0.416  = AU] 


across  flame  before 
Interaction 


across  the  flame  front  after 
Interaction 


Hence  the  aolutuion  is  OK 
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Interaction  of  the  Flame  front  with  transmitted  shock  in  the  balancing  chamber:— 


Figure  (13)  shows  the  interaction  involving  the  flame  front  and  a transmitted  shock  of 
strength  given  by  Ms  = 1.258.  The  matching  condition  across  the  flame  front  are  the 
pressure  and  flow  velocity. 


/ 

/ 


Fig.  13.  interaction  of  the  shock  Wave  with  a Flame  Front: 


Guess  Ajq  2.21 

Then  ~~  = 1.0549 
An 


Now  from  the  shock  tables,  we  have 
Ms  = 1.177 


AS  - 0.004 


I A U I 

A 


= 0.273 


etc. 


0 
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The  expansion  wave  resulting  from  this  interaction  is  of^onsideraNe  significance. 
Thus  it  is  interesting  to  note  that  flow  velocities  in  regions  (z^  and  (3j)  across  the 
expansion  wave  are  of  opposite  signs.  Hence  there  is  a change  in  flow  direction  across 
the  expansion  wave.  The  interface  also  changes  direction  as  a result  of  this  expansion 
wave  and  accelerates  towards  the  main  combustion  chamber  and  finally  reaches  the 
flame  front  at  point  A when  combustion  stops.  The  interface  thus  travels  towards  the 
main  chamber  from  the  balancing  chamber  as  a result  of  wave  interactions. 

Thus  the  slowing  down  of  the  flame  front  and  the  movement  of  the  interface 
towards  the  main  chamber  is  of  extreme  importance  because  of  the  availability  of  oxygen 
towards  the  end  of  combustion  and  the  increased  time  available  for  completion  of 
combustion  process  in  the  NAH6E  engine. 


